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Summary. — Recent astronomical observations of SNIa, CMB, as well as BAO 
in the Sloan Digital Sky Survey, suggest that the current Universe has entered a 
stage of an accelerated expansion with the transition redshift at z ~ 0.5. While the 
simplest candidates to explain this fact is cosmological constant/vacuum energy, 
there exists a serious problem of coincidence. In theoretical cosmology we can find 
many possible approaches alleviating this problem by applying new physics or other 
conception of dark energy. We consider state of art candidates for the description of 
accelerating Universe in the framework of the Bayesian model selection. We point 
out advantages as well as troubles of this approach. We find that the combination 
of four data bases gives a stringent posterior probability of the ACDM model which 
is 74%. This fact is a quantitative exemplification of a turmoil in modern cosmology 
over the A problem. 

PACS 95.36. +x - Dark energy. 

PACS 98. 80. Es - Observational cosmology. 



1. Introduction 

Recent observations of SNIa stars indicate that the Universe is in an accelerating phase 
of expansion. If we assume that the FRW model with a source in the form of a perfect 
fluid well describes present evolution of the Universe, then there is only way to explain the 
observational fact that the Universe is accelerating, i. e. postulating some gravity source of 
unknown nature which violates the strong energy condition p + 3p > 0. The most natural 
and simplest candidate for such a type of dark energy is the vacuum energy represented 
by the cosmological constant A in the ACDM model. But we have still a problem with 
this model: predicted value of the vacuum energy density is larger than what we observe 
by a factor of the order 10 120 (the so called fine tuning problem) and energy density 
parameters of both dark matter and dark energy are comparable at the present epoch 
(the so called coincidence problem) . Due to that we are still looking for other alternative 
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Table I. - The Hubble function for cosmological models with dark energy 



model 



W* 



ACDM model 

Model with generalized 

Chaplygin gas 

Model with phantom 

dark energy 

Model with dynamical 

E.Q.S 

Quintessence model 



M 1 + z f + C 1 - n m,()) 
■ n m , (l + zf + (1 - n m , )[^s + (1 - A S )(1 + Z )3(1+«)]TTS 

n m , (i + zf + (i - n m , )(i + zf^+ w x) 

n m , (i + + (i - n m , )(i + s )3C«o+»i+i) exp [- 3 ^f] 

n m , (i + zf + (i - n m , )(i + 2)3(i+«. (i+*)- a ) 



explanations of accelerating expansion phenomenon. Apart from the dark energy models 
we have a second group of models where the modification or generalisation of the FRW 
equations are postulated. We gather the most popular models belonging to the sets 
deseribed before and we try to answer the question if they are better than the ACDM 
model in the light of available data sets using the Bayesian model comparison method. 



2. Cosmological models of the accelerating Universe 

We consider five models belonging to the group of models with dark energy and gather 
them together with Hubble funetions, with the assumption that the Universe is spatially 
flat in Table HJ Here we include the ACDM model [T], the phantom dark energy model 
[2j [3], where the coefHcient of the equation of state is a parameter (w — const) with 
the prior assumption that it is less than —1, a model with dynamical equation of state 
[4l [5] where the coefHcient of the equation of state is a linear function of a: w(a) — wq + 
u>i(l — a). We also consider the quintessence model where dark energy is represented by a 
minimally coupled scalar field 0(7] . We choose the power-law parameterized quintessence 
model, where the evolution of the dark energy density is given by px — px.o a ^ 3< ' 1+Wx< " a ^ 
and wx (a) is the mean of a coefHcient of the equation of state in the logarithmic scale 
faetor w)x (a) = / wx{a)d{lna)/ J d(lna) with parametrization proposed by [8j : wx(a) = 
woa a . Finally we include model with the generalized Chaplygin gas [9j [TOl [Tl] , i. e. with 
perfect fluid which equation of state has the following form px — —A/p x , where A > 
and a is a constant. 

We also consider five models which belong to the the group of models with the modified 
or generalized FRW equation (gathered together with the Hubble funetions in Table lH|) . 
Here we include examples of so called brane models, which postulated that the observer 
is embedded on the brane in a larger space in which gravity can propagate: the Dvali- 
Gabadadze-Porrati model (DGP) [H], the Sahni-Shtanov brane 1 model [HC3]. Thcrc 
is also the Cardassian model where the modification of the Friedmann first integral is 
postulated: ZH 2 = p + Bp n , and S is a constant |15j . We also consider the bouncing 
model arising in the context of loop quantum gravity (the BACDM model) [UI [17] and 
the model with energy transfer between the dark matter and dark energy seetors |18) . 
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model 



H' (z) 



DGP 

BACDM 

interacting 
with A 

Cardassian 

n r , = icr 4 

Sahni-Shtanov 
branc I 



y/ fV,o(l + z) 3 + Q r c,0 + y/Clrcfi 
_ (1-O m ,o) 2 



n m ,o(i + zf- - n n>0 (i + z) n + 1 - n m , + n„ i0 
fim,o(i + z ) 3 + ^mt,o(i + z) n + 1 - n m , - n Jnt ,o 



n r ,o(i + ^) 4 + n m ,o(i + 2:) 4 



^ + (l + z) -4 + 4n( l-y-..o ^ (2) 



= n m , (i + z) 3 + n CT , + 2c ; ,o - 2 v /n^p(«) 
p(z) = ^/n mfi (i + Z) 3 + 

,0 + ^1,0 + ^Af>,0 



3. Bayesian model comparison framework and application to cosmological 
models 

In the Bayes theory the best model from the set under consideration is the one which 
has the largest value of the posterior probability [19] (the probability in the light of data 
D) which has the following form 

d, PmD) = ffiwp . 

P(M) is the prior probability for the model considered, the value of which depends on 
our previous knowledge (i. e. without information coming from the data D) about the 
model. If we have no foundation to favour one model over another one from the set 
of models under consideration we usually assume the same value of the prior quantity 
for all models, i. e. P(M) — l/K, where K is the number of models. P{D) is the 
normalization constant. Requiring that a sum of posterior probabilities for all models 
from the considered set is equal to one we can obtain the form for this quantity: 

K K 

(2) ^P(Mi\D) = l -^P(D)=J2P(D\M i )P(M i ), 



where i indexed the model under consideration. It should be pointed out that conclusions 
based on the posterior probability depend on the set of models under consideration and 
could change if the set is changed. 

P(D\M), the marginal likelihood (also called the evidence), is the most important 
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quantity in the Bayesian framework of model comparison and has the following form 



(3) 



P(D\M) = / L(6\D,M)P(0\M)d0 = E, 



where L(9\D 1 M) is the likelihood of the model M, 9 is the vector of the model parameters 
and P(6\M) is the prior probability for the model parameters. Here we use an asymptotic 
approximation to the minus twice logarithm of the evidence derived by Schwarz [20], so 
called BIC quantity which has the following form: 



(4) 



BIC 



-21n£ + dln7V, 



where C is the maximum of the likelihood function, d is the number of model parameters 
and N is the number of data used in the analysis. This approximation is derived with the 
assumption that observations {D = {xi}, i = 1, • • • , N) used in analysis are iid (indepen- 
dent and identically distributed) and comes from a linear exponential distribution defined 



in the following way f(xi\9) 
td 



exp 



where wi, ■ ■ ■ , w s, b are 

functions of only 6 S R d , t\,. .. ,ts are function of only Xi and S = d. The prior for 
model parameters must be not equal to in the point where the likelihood function reach 
maximum and the number of data must be large enough with respect to the number of 
model parameters (as a consequence of asymptotic approximation) . This derivation was 
generalized by [21] where it is not required to assume any specific form for the likelihood 
function but is only necessary that the likelihood function satisfies some non-restrictive 
regularity conditions. Moreover the data do not need to be iid. This derivation requires 
to assume that a prior for model parameters is not equal to zero in the neighborhood 
of the point where the likelihood function under a given model reaches a maximum and 
that it is bound in the whole parameter space under consideration. The number of the 
data must be also large with respect to the number of model parameters. 

We apply this method to compare the cosmological models described before. We 
make the comparison among the models with dark energy (gathered in Table[I|, among 
the models with modified theory of gravity (gathered in Table|TT]) as well as among all 
models described before. We use the information coming from the SNIa, CMB, BAO 
and H(z) data. Let us describe those data sets. 

We use Ni = 192 SNIa data [22] [23l [24] . In this case the likelihood function has the 
following form 



(5) 



L sn oc exp 



g (/4 h cor 



obs\2 



where Oi is known, /i° bs = mi — M (mj-apparent magnitude, M-absolute magnitude 
of SNIa), /4 heor = 51og 10 D u + M, M = -51og 10 ff + 25 and D Li = H d Lil where 
d^i is the luminosity distance, which with assumption that k = is given by d^ — 
(l + Zl )cf z '(dz>/H(z>)). 

We also include information coming from the CMB data. Here the likelihood function 
has the following form 



(6) 



Lr oc exp 



(K 



thcor 



obs\2" 



R ohs ) 



24 
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Table III. - Posterior probabilities for models from Tablet and for models from Table UT\ 



model 


prior 


posterior 


model 


prior 


posterior 


1 


0.20 


0.84 


6 


0.20 


0.07 


2 


0.20 


0.02 


7 


0.20 


0.03 


3 


0.20 


0.06 


8 


0.20 


0.13 


4 


0.20 


0.04 


9 


0.20 


0.74 


5 


0.20 


0.04 


10 


0.20 


0.03 



where R is the shift parameter, i? theor = y/H^ f* d ° c (H /H(z))dz, and R ohs = 1.70 ± 
0.03 for 2dec = 1089 [2S1I2B]. 

As third observational data we use the measurement of the baryon acoustic oscillations 
(BAO) from the SDSS luminous red galaxies [27] ■ In this case the likelihood function 
has the following form 



(7) 



La oc exp 



(A 



thcor 



2o\ 



where A thcm = y/U^ (H(z A )/H )- 1/3 [(l/z A ) f* A (H /H(z))dz] 2/3 and A ohs = 0.469± 
0.017 for z A = 0.35. 

Finally we use the observational H(z) data (N 2 — 9) from [28] (see also [29] [30] and 
references therein). These data are based on the differential ages (dt/dz) of the passively 
evolving galaxies which allow to estimate the relation H{z) = d/a = —1/(1 + z)(dz/dt). 
Here the likelihood function has the following form 



Lu oc exp 



' N 2 

E 

i=l 



where H (z) is the Hubble function, Hi, Zi are observational data. 

The final likelihood function is given by L = LsnLrLaLh, with the number of data 
points = 192 + 1 + 1 + 9. We assume equal prior probabilities for all models from 
the set considered. Results, i. e. the values of the posterior probabilities, for the set of 
models with dark energy and set of models with modified theory of gravity are gathcred 
in Table [m] whilc for the set of all models in Table [IV] 

We can conclude that the ACDM model is the best one in the set of model with dark 
energy as well as the best one among all considered model. If we compare only models 
with modified theory of gravity the best one from them is the Cardassian model. 



4. Conclusion 



Our results of Bayesian method of model selection indicate that the ACDM is the best 
one. On the other hand this model has the status of an effective theory which gives us a 
good dcscription, however without deep understanding. The ACDM model is based on 
the general relativity theory which is believed to be superposed by a more fundamental 
theory, namely quantum gravity theory. Therefore, we still search for a quantum based 
model of the Universe [321 , 33J . So far the ACDM model is some effective theory of the 
Universe the best we have. In our opinion the ACDM model is expected to be emergent 
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Table IV. - Posterior probabilities for models from Tables\]\ and [771 



model 


prior 


posterior 


1 


0.10 


0.74 


2 


0.10 


0.02 


3 


0.10 


0.05 


4 


0.10 


0.04 


5 


0.10 


0.03 


6 


0.10 


0.01 


7 


0.10 


0.005 


8 


0.10 


0.01 


9 


0.10 


0.09 


10 


0.10 


0.005 



from thc future quantum gravity theory of the Universe and finding it is thc challcnge of 
the 21st century. 
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